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We have performed spectroscopic measurements of a superconducting qubit dispersively coupled 
to a nonlinear resonator driven by a pump microwave field. Measurements of the qubit frequency 
shift provide a sensitive probe of the intracavity field, yielding a precise characterization of the 
resonator nonlinearity. The qubit linewidth has a complex dependence on the pump frequency and 
amplitude, which is correlated with the gain of the nonlinear resonator operated as a small-signal 
amplifier. The corresponding dephasing rate is found to be close to the quantum limit in the low-gain 
limit of the amplifier. 



PACS numbers: 85.25.Cp, 74.78.Na, 03.67.Lx 

In quantum mechanics, any measurement necessar- 
ily induces decoherence in the variable conjugate to the 
one being measured. This principle is quantitatively ex- 
pressed by the inequality r^, m > r moas /2 between the 
system's measurement-induced dephasing rate r^, m and 
the measurement rate r mcas , stating that the most effi- 
cient detector can only measure as fast as it dephases [1]. 
Determining how far a specific detector departs from this 
quantum limit is a fundamental issue. The measure- 
ment of a two-level atom by a cavity to which it is non- 
resonantly (dispersively) coupled has been studied in de- 
tail theoretically and experimentally, first in cavity quan- 
tum electrodynamics (QED) with Rydberg atoms [3], 
then in circuit QED with superconducting qubits [4- 
6]. In this situation, information about the atom state 
i (ground state g or excited state e) is encoded in the 
complex dimensionless amplitude of the qubit-dependent 
intracavity field a; (the pointer states [2]) and is quanti- 
fied by the distinguishability D = \a e — a s \. The back- 
action of this dispersive readout consists in a shift of the 
atom frequency called the ac-Stark shift and in a gradual 
dephasing at a rate r<(, m = kD 2 /2 (k is the resonator 
damping rate) reaching the quantum limit r moas /2 [6]. 

Introducing a nonlinearity into the measuring cavity 
makes it possible to turn it into an active device, which 
can considerably enhance the signal-to-noise ratio, a key 
requirement for high-fidelity measurement of supercon- 
ducting qubits [7-9]. In appropriate parameter ranges, 
pumped nonlinear resonators indeed behave as noiseless 
parametric amplifiers for small incoming signals [10], gen- 
erate squeezing [11], and display bistability [12-14]. An 
important fundamental open question is then to deter- 
mine in what respect these novel physical effects modify 
the backaction of the measuring cavity onto the atom, 
and to investigate whether the measurement is still at 
the quantum- limit or not. 



To address this issue, we measure the spectrum of a su- 
perconducting transmon qubit [15] dispersively coupled 
to a nonlinear resonator pumped by a microwave tone. 
Both the ac-Stark shift and the measurement-induced de- 
phasing show new features compared to the case of a 
passive linear resonator. The ac-Stark shift provides a 
sensitive probe of the intrarcsonator mean photon num- 
ber n and allows us to characterize the resonator non- 
linearity with high precision. The measurement-induced 
dephasing rate is found to have a completely different de- 
pendence on the pump power than in a linear resonator: 
instead of growing linearly with n, T^, m peaks where the 
parametric amplifier gain is maximum and decreases to- 
wards zero at large n. This unexpected behavior is well 
described by the simple relation — kD 2 /2, D being 
now the distinguishability between the pointer states of 
the nonlinear resonator. We finally discuss how close this 
backaction is to the quantum limit. 

In our experiment [see Fig. la)], a transmon qubit is ca- 
pacitively coupled with a strength g/2n = 42.4 MHz to a 
A/2 coplanar resonator of frequency cj /2tt — 6.4535 GHz 
and quality factor Q = 670 (damping rate k/2tt = 
9.6 MHz). The resonator is made nonlinear by insert- 
ing a Josephson junction of critical current Io = 750 nA 
in its center [9]. We model [16, 17] this nonlinear dis- 
tributed resonator by one single mode with Hamiltonian 
H/h = LO ^a+(K/2)(ai) 2 a 2 + (K'/3)(ai) 3 a 3 containing 
a Kerr nonlinearity with constants K and K' = 2-lO~ 3 K 
that can be explicitly derived from the circuit parameters 
(see Supplementary Material) . The classical steady-state 
response a of such a Kerr resonator (KR, also called Duff- 
ing oscillator) to a pump drive of frequency w p and input 
power P p is therefore highly nonlinear and given by [18] 

i (n^a + K\a\ 2 a + K'\a\ 4 aj + = -ie p , (1) 

where = 2Q(1 — uj p /uj ) is the reduced detuning of the 
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FIG. 1. (color online) (a) Experimental setup: a transmon 
qubit is strongly coupled to a coplanar resonator made non- 
linear with a Josephson junction. The sample is cooled to 
20 mK and is driven through an attenuator A by three mi- 
crowave sources. Source p is used to establish a pump field in 
the resonator, source q for qubit spectroscopy, and source r as 
a JBA readout: its signal at uj t is reflected from the resonator 
and routed through circulators to a cryogenic amplifier, a de- 
modulator, and a digitizer, which yields the JBA switching 
probability p s and thus the probability of the qubit excited 
state e. (b) Frequency dependence of the maximum current / 
in the resonator calculated for the sample parameters and for 
reduced powers P P /P C = 0.1,0.635, 1,3.5, and 16.35 (bottom 
to top), (c) Stability diagram of the resonator in the f2-P p 
plane. The solid lines indicate the highest PA gain below 
Q c and the power at which the resonator bifurcates from the 
low- (L) to the high- (H) amplitude state, respectively. The 
hatched area is the bistability region. Vertical dotted lines 
correspond the two datasets of Figs. 2-3. 



pump and e p = ^KP p /Ahaj p its reduced amplitude (A is 
the total attenuation of the input line). The maximum 
amplitude of the oscillating current / = h / ir ZoUJo\a\ 
in the resonator is shown for our sample parameters 
in Fig. lb). At low drive amplitude, the response is a 
Lorentzian around £1 = as for a linear resonator. For 
stronger drive, the resonance frequency shifts downwards 
and shows a sharpened response in a frequency window 
in which the resonator behaves as a parametric amplifier 
(PA) [10], until the slope becomes infinite for a critical 
power P c at f2 c = V3. For P p > P c and f2 > Q c , two 
stable solutions of different oscillation amplitude can co- 
exist, L (low) and H (high). In this bistable regime called 
bifurcation amplification (BA), the transition from L to 
H occurs abruptly when ramping up the pump power 
at the bifurcation threshold P+(fi). Fig. lc) summarizes 
these properties. 
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FIG. 2. (color online) 2D plots of the switching probability 
p s (uj q ,P p ) for fi g /fi c = 3.1 and 0.7 (left and right panels). A 
few qubit lines are shown in overlay for negligible field am- 
plitude in the resonator (top left), near the switching point 
at P p = 1.0 dBm for fi g /n c = 3.1, and at P p = -12.4 dBm 
and —8.6 dBm for Q s /Sl c — 0.7. The horizontal dashed line 
indicates the qubit frequency at zero cavity field. Lorentzian 
fits of the qubit lines (see example at —8.6 dBm) yield the 
ac-stark shift Aoj q and the FWHM linewidth w. Inset: mi- 
crowave pulse sequence used. 

The qubit-resonator detuning is fixed at A/2n = 
732 MHz 3> g/2n, so that their interaction is well de- 
scribed by the dispersive Hamiltonian H mt = h(\a z + 
s)a^a with x/2tt = -0.8 MHz and s/2tt = 1.7 MHz [19]. 
In this dispersive regime, the resonator frequency takes 
the qubit-state dependent value Wj = luq + s ± \ for 
i = g,e. Consequently, the reduced pump detuning 
Hi as well as the bifurcation threshold P+(f2;) now de- 
pend on the qubit state. This allows us to readout 
the qubit by sending microwave pulses with a frequency 
w r /27r = 6.439 GHz (f2 g /r2 c = 2) and power P r between 
P + (ri e ) and P+(i7g) [9]. After reflection on the cavity, 
the phase of the readout pulse is measured by homo- 
dyne detection yielding the final oscillator state. By re- 
peating this sequence, the Josephson Bifurcation Ampli- 
fier (JBA) switching probability p s , and thus the qubit 
excited-state probability, is determined. 

To study the backaction of the pumped KR, we per- 
form qubit spectroscopy with a microwave pulse of vary- 
ing frequency uj q and fixed power P q while the cavity is 
driven by a microwave pump pulse of varying frequency 
ujp and power P p . The experimental sequence ends with 
a qubit readout pulse (see inset of Fig. 2). Note that 
the pump pulse starts long before the spectroscopy pulse 
so that the intracavity field has reached its stationary 
state. Moreover, the readout pulse is applied 200 ns af- 
ter switching off both other pulses, a time long enough 
to let the intraresonator field relax before readout, but 
shorter than the qubit relaxation time T\ = 700 ns. Fig. 2 
shows the resulting qubit spectrum as a function of P p 
at two pump frequencies above and below f2 c . We ob- 
serve the ac-Stark shift of the qubit resonance frequency 
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towards lower values, as well as the broadening of the 
qubit line. This is akin to the behavior observed with 
a linear resonator, but with a very different dependence 
on P p [5]. First, an abrupt discontinuity in the ac-Stark 
shift at J7 g /£l c = 3.f clearly indicates the sudden in- 
crease in the intracavity average photon number n as 
the resonator switches from L to H. In this region, two 
spectroscopic peaks are observed at a given pump power. 
Second, the linewidth narrows down at large n, in strong 
contrast with the linear resonator case where it increases 
linearly [5]. 

Our measurements of the qubit frequency shift Aw q 
versus P p allow us to quantitatively determine n, us- 
ing a model for the ac-Stark shift that includes the 
next order correction to the linear ac-Stark shift for- 
mula n — Aoj q / (2x) [20] and that takes into account 
the first 5 levels of the transmon [15, 19] (see Fig. 3a). 
We then fit the resulting experimental n(P p ) curves with 
the n values calculated from the sole dynamics of the 
resonator (with the qubit in g) given by the square mod- 
ulus of the solutions of Eq. (1) using K and A as the 
only fitting parameters. The agreement is excellent over 
the whole (ft, P p ) range for K/2ir — —625 ± 15 kHz and 
^4 = 110.8±0.2 dB, which is consistent with the design 
value K/2tt — —750 ± 250 kHz and with independent 
measurements of the line attenuation A = 111 ± 2 dB. 
This demonstrates that measuring the ac-Stark shift of 
a qubit is a sensitive method for probing the field in- 
side a nonlinear resonator and for characterizing its Kerr 
nonlinearity [10]. 

We show in Fig. 3b) the qubit total dephasing rate 
r 2 = 7rw, with w the full width at half maximum 
(FWHM) of a Lorentzian fit to the experimental data. 
In addition to the field- induced dephasing T^m, r 2 also 
includes a constant contribution r 2 due to all other de- 
phasing processes (mainly energy relaxation, flux noise 
and radiative broadening). For the sake of comparison, 
we also show the independently measured gain G of the 
KR operated as a parametric amplifier on a small addi- 
tional signal frequency-shifted by 100 kHz from uj p /2tt. 
Below and above fi c , r 2 (P p ) peaks where G is maximum 
and near the bifurcation threshold respectively, and tends 
towards r 2 at large n. 

To account for this behavior, we have calculated r 2 
by numerical integration of the master equation of the 
multi-level transmon coupled to the KR with a Jaynes- 
Cummings Hamiltonian for the sample parameters and 
damping rates determined independently [see full sym- 
bols in Fig. 3b)]. The agreement with the experimental 
data is excellent with no adjustable parameter, in con- 
trast with experiments with a flux-qubit in which an un- 
explained shortening of the qubit relaxation time is ob- 
served when the KR is driven above bifurcation [21, 22]. 
We also derive an analytical expression for r^ m starting 
from the master equation in the dispersive limit, using 
methods similar to those of [20] for the linear case. This 
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FIG. 3. (color online) (a) Experimental (open symbols) av- 
erage photon number n and corresponding fits (lines - see 
text) as a function of P p for the same datasets n g /fi c = 3.1 
(circles) and n g /Sl c = 0.7 (triangles) as in Fig. 2. Experi- 
mental points are obtained by converting the measured Auu q 
according to the ac-Stark shift model of the inset (solid line 
- see text ; the dashed line shows the linear approximation), 
(b) Qubit dephasing rate T2 measured for the same datasets 
(open symbols) and calculated either by numerical integra- 
tion of the system master equation (solid symbols) or using 
Eq. (2) (solid lines). The horizontal dashed-dotted line in- 
dicates the intrinsic dephasing at zero field. The measured 
parametric power gain (dashed line) is also shown for com- 
parison for Q s /Sl c = 0.7. (c) Complex cavity field amplitude 
a g (solid and dashed lines) calculated from Eq. (1) for in- 
creasing pump powers and for the same datasets. Both a g 
(squares) and a e (triangles) are also shown at six points la- 
belled A-F. Their separation (segments) is to be compared to 
the uncertainty disk of a coherent state (open circles or disks 
shown at point C and at points A-F in panel b). 



is done by linearizing the quantum fluctuations of the 
field around the two classical solutions a\ of Eq. (1) that 
correspond to the two qubit states, assuming that the 
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separation between these pointer states is small and ne- 
glecting squeezing and transients (T^ <C k). The details 
of these calculations, which go beyond a simple linear re- 
sponse approximation, will be presented elsewhere [23]. 
In this way, we obtain 

IV = f|a -a g | 2 = ^ 2 . (2) 

We thus recover the same link as for a linear resonator 
between decoherence rate and distinguishability. Here 
however, due to the nonlincarity, D and therefore 
are no longer proportional to n [5] and on the contrary 
tend towards zero at large n. Since in a pumped KR the 
susceptibility \da/dfl\ varies similarly to the small-signal 
gain G [24, 25], Eq. (2) also explains the correlation ob- 
served between and G. The prediction of Eq. (2) 
calculated without any adjustable parameter as well as 
the corresponding pointer states are shown in Fig. 3b) 
and c). The theory reproduces well the non-trivial de- 
pendence of Y2 on the pump power over a large parame- 
ter range, and in particular the dephasing peak observed 
above and below f2 c . The direct correlation between T 2 
and D is clear. The agreement is quantitative in the 
validity range of our approximation D<1, which is sat- 
isfied by our data when T 2 is below ~ 10 MHz. 

In order to verify whether the pumped KR backaction 
is at the quantum limit, needs to be compared to 
the rate r meas at which information about the atom state 
leaks out of the resonator. If the two pointer states are 
coherent states, r mcas = kD 2 [6] and our results therefore 
establish that the quantum limit is reached. According 
to the quantum theory of the pumped KR [25, 26], this 
is however true only in the limit of small parametric gain 
G <~ 1 [corresponding to P p significantly different from 
P+(f2)]. Indeed, when G ^> 1, the intracavity field is 
expected to show enhanced phase-dependent fluctuations 
and the field reflected on the cavity to show some degree 
of squeezing. The measurement rate then needs to be re- 
evaluated. Note that recent theoretical work [25] showed 
that if the linear response theory is valid and if the gain 
is large, the backaction of the pumped KR misses the 
quantum limit by a large factor of order G. In contrast 
our experimental parameters, which are typical for circuit 
QED, require going beyond linear response but involve 
limited parametric gain. It is also important to note that 
our work docs not address the situation encountered in 
the BA regime close to the bifurcation threshold P + (fl), 
in which a gje are the two different oscillator states L, H. 
As a result, further work is needed to decide whether the 
readout of a qubit by a JBA is quantum- limited or not. 

In summary, we have measured the quantum backac- 
tion of a pumped nonlinear resonator on a dispersively 
coupled qubit. We observe a non-trivial dependence 
of the measurement-induced dephasing on the pump 



power, which we link to the distinguishability between 
the pointer states of the nonlinear resonator and to the 
gain of the resonator operated as a parametric ampli- 
fier. In the small gain limit, we find that the backaction 
is quantum-limited; further theoretical work is needed 
to describe the high gain regions. More generally, our 
results demonstrate that circuit QED is an ideal play- 
ground to study the interplay between strong coupling 
and nonlinear effects such as bistability, squeezing or 
parametric amplification. 
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In this Supplementary Material, we obtain analyti- 
cal expressions for the Kerr constants appearing in the 
Hamiltonian of a nonlinear resonator. It consists of a A/2 
coplanar waveguide resonator with an embedded Joseph- 
son junction of critical current Jo embedded in the mid- 
dle conductor. We first discuss why it is appropriate 
to model the distributed nonlinear resonator by a single 
mode with Kerr nonlinearities. Indeed, in doing so we ne- 
glect the higher frequency modes of the resonator. This 
is justified in our case because all the resonator pumping 
is performed very near the Kerr resonator frequency ujq, 
and because the resonance frequency of the other modes 
of our resonator is far from being a multiple of ljq as dis- 
cussed in [1], allowing us to neglect cross-Kerr coupling 
between various modes [2] . 

To derive the Kerr constants, the first step is to map 
the distributed nonlinear resonator onto an equivalent 
series combination of a lumped element inductance L e , 
capacitance C e (which includes the junction capacitance) 
and Josephson junction of critical current Ig, as shown 
in Fig. 1. Note that this mapping is only possible here 
because we can neglect cross-Kerr terms between the var- 
ious modes of the resonator sin In the limit where the 
Josephson inductance Lj is completely negligible com- 
pared to the resonator inductance, one can show that 
L c = 7tZo/2oji and C c = 2/wZqUJi, where Z 50 ft is 
the resonator characteristic impedance and uj\ the res- 
onance frequency in absence of the junction. When Lj 
is not negligible, one has to adjust numerically L c and 
C e for the impedance of the equivalent circuit to fit the 
exact distributed resonator impedance. This is the ap- 
proach that has been used in the Letter. 

Using the notation introduced in Fig. 1, we obtain the 
equivalent circuit Hamiltonian 



H 



2L r 



Ej cos 



<r 

2C e 



(1) 



Since the current I flowing through the inductance and 
the junction is the same, we also have 



0i 

I = — = J sm 



9o 



(2) 



yielding an implicit relation (\>\ — g(<f>) between the two 
phases. Eliminating <f>\, the Hamiltonian thus takes the 
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FIG. 1. Equivalence between the distributed nonlinear res- 
onator and a series combination of an equivalent inductance 
Le, capacitance C e and Josephson junction of critical current 
lo- 



form 
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By expanding g(ifi) in powers of <fr, we obtain the non- 
linear resonator Hamiltonian to any order of Josephson 
junction nonlinearity. For instance, to fourth order we 
find 



H 



4? 

2Lt 



q 

2C C 



2A P 



(4) 



where L t = Lj + L c is the total inductance and p = 
Lj/Lt. 

This Hamiltonian can be writen in terms of the 
creation and annihilation a operators with <f> — 
iyfhZ c /2(a — a)) and q — y // h/2Z c (a + a 1 ) with Z c = 
\J Lt/Cc- Once expanded, the nonlinear term cj} 4 yields 
products of creation and annihilation operators to vari- 
ous powers. Using the rotating-wave approximation, we 
keep only those with equal number of annihilation and 
creation operators yielding 



K 

H = hu)Q(ra + h—(a)) 2 a 2 , 



where 



K = - 



TTp 3 0J Z e 

Rk 



(5) 



(6) 



is the Kerr constant, cu = l/\/L t C c and Rk — h/e 2 . 
Here and below, we have dropped the small correction to 
Wo coming from K. 
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In the same way, the Hamiltonian can be expanded 
to the next order of junction nonlinearity yielding the 
next-order Kerr constant K' 



with 



K K' 
H = huoa'a + h—{a}) 2 a 2 + h—(a r ) 



2 K 2 

K' = ~ (10p-9) 
3p w 



(7) 



(8) 



relation K' w 0.002.K" that we consider fixed throughout 
the Letter and in particular in the fit of the ac-Stark shift 
data. 
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For our sample's parameters, this yields an approximate 



